Density-relaxation part of the self energy 

In a recent Letter, Ogiit, Chelikowsky and Louie Jl]] 
presented an important series of calculations of the ef- 
fect of quantum confinement on optical gaps in large 
hydrogen-passivated spherical silicon clusters, by calcu- 
lating the quasiparticle energies for addition of an elec- 
tron and of a hole separately, and then the excitonic bind- 
ing energy. The quasiparticle energies were calculated 
by what might be termed a ALDA approach: within the 
local-density approximation (LDA) , the ground-state to- 
tal energies of the n-, n — 1- and n + 1-electron systems 
(where n is the number of electrons in the neutral clus- 
ter) were calculated, and then the quasiparticle gap was 
estimated using 

ef = E n+1 + E n _ x - 2E n . (1) 

The authors suggest that this expression would be ex- 
pected to approach the experimental quasiparticle en- 
ergy gap of bulk silicon (1.2 eV) in the large-cluster 
limit. They presented a numerical fit of the correction 
e g P ~ e g P LDA' wn ich they stated approached the bulk value 
of 0.68 eV like where £g P LDA is the LDA Kohn- 

Sham eigenvalue gap and d is the cluster diameter. How- 
ever, it is known that in the bulk limit Eq. (|l|) (in the 
LDA) simply yields the LDA energy gap: the correction 
is zero. This is because the LDA exchange-correlation 
energy is an analytic functional of density: the fact that 
the change in electron density on adding (or subtract- 
ing) a single electron is of order 1 /n allows the changes 
in the Kohn-Sham eigenvalues and the other ingredients 
of the energy to be evaluated using perturbation theory, 
and after a substantial cancellation between terms the 
stated result is obtained. (The same formula yields the 
correct gap in exact Kohn-Sham DFT, but this reflects 
a non-analytic discontinuity in the exchange-correlation 
potential between the n- and n+ 1-electron systems |Q.) 

In physical terms, the ALDA approach includes the 
electrostatic effect of the relaxation of the charge density 
when an electron is added or subtracted, and the cor- 
responding relaxation in the LDA exchange-correlation 
potential. In the large-cluster limit, both these effects 
go to zero, and the non-zero band-gap correction may 
be calculated using many-body perturbation theory in a 
suitable approximation (e.g. I^Q), where the correction 
to the LDA band gap arises from the differing effects of 
the non-local self-energy on the states concerned || . 

Furthermore, there is no reason to suppose that this 
term in the self-energy correction that is excluded in 
the ALDA approach is negligible in the clusters stud- 
ied. Therefore it is likely that the quasiparticle and op- 
tical gaps given in Ref. Q should be increased by very 
approximately 0.68 eV, where the error bar in this esti- 
mated correction is smallest for the largest clusters. Of 



course, this additional correction is of lower relative im- 
portance for the smaller clusters. 

To confirm our theoretical analysis, we have reanalyzed 
the data for the ALDA gap correction from Ref. [0, 
shown here in Fig. 1 as a function of 1/d. The 
dashed curve shows the best (least-squares) fit of the 
form 0.68 eV + Ad~ p , as in Ref. [0, obtained by us 
with p =1.40 (similar to the 1.5 given in Ref. 0). The 
solid curve shows the best fit obtained if the constraint 
that the limit as d — > oo should be 0.68 eV is removed 
(as it should be): K + Ad~P with K =(0.12±0.24) eV, 
p = 0.92 ±0.14. The value of K is indeed consistent 
with zero, and inconsistent with 0.68 eV. The second fit 
is more than twice as good as the first, as measured by 
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FIGURE CAPTIONS 

Fig.l The band-gap correction from Ref. [Q, plotted 
against the inverse cluster diameter. The best fit 
(solid curve) correctly tends to a value consistent 
with zero. 
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